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PERFORMANCE AND CONTROL STUDY OF A LOW-PRESSURE-RATIO 
TURBOJET ENGINE FOR A DRONE AIRCRAFT 


by Kurt Seldner, Lucille C. Geyser, Harold Gold, 

Darrel Walker*, and Gary Burgner* 

Lewis Research Center 
SUMMARY 

The results of analog and digital computer studies of a low pressure -ratio turbojet 
engine system for use in a drone vehicle are presented. The turbojet engine consists 
of a four-stage axial compressor, single-stage turbine, and a fixed area exhaust nozzle. 
Three simplified fuel schedules and a generalized parameter fuel control for the engine 
system are presented and evaluated. The evaluation is based on the performance of 
each schedule or control during engine acceleration from a windmill start at Mach 0. 8 
and 6100 meters (20 000 ft) to 100 percent corrected speed. It was found that, because 
of the higher acceleration margin permitted by the control, the generalized parameter 
control exhibited the best dynamic performance. 


INTRODUCTION 

The low-pressure-ratio turbojet engine is a suitable propulsion system for many 
drone vehicles. The short flight range of most drone missions places minimal require- 
ments on engine fuel consumption. The basic engine simplicity that accompanies the ** 
low pressure ratio makes possible low-cost engine construction. Low engine cost is, 
in turn, a requirement for the economic feasibility of many drone applications. 

This report presents the results of a performance and control study of a four to one 
compression ratio engine operating in a drone vehicle cruising at Mach 0. 8 at 6100 me- 
ters (20 000 ft) and having a thrust requirement of 1557 newtons (350 lbf). The vehicle 
was assumed to be started and air launched at the cruise conditions. The design analy- 
sis presented was performed on analog and digital computers and is based on the meth- 
ods described in references 1 and 2. 

*Naval Weapons Center, China Lake, California. 



This study was made under the Lewis Research Center low-cost engine study pro- 
gram and was performed in cooperation with the Naval Weapons Center, U. S. N. 

Design Point Characteristics 

In the initial phases of the study, a preliminary cycle analysis based on the com- 
puter program of reference 3 was performed on the engine system. This study was used 
t<j establish a suitable engine design operating point. Results obtained from this analy- 
sis are plotted as variations in specific thrust against specific fuel consumption for sev- 
eral compressor pressure ratios and a range of turbine inlet temperatures in figure 1. 

A final selection of a pressure ratio of 4 and a turbine inlet temperature of 1078 K 
(1941° R) was based on the simplicity of the design needed to accommodate system oper- 
ation at this condition. 

Table I presents the actual design data and computed system parameters for various 
operating conditions. The compressor efficiency and rotational speed are the design 
values for a four-stage axial-flow transonic compressor. The turbine efficiency is the 
design value for a single-stage axial-flow machine. The combustor outlet pressure cor- 
responds to an annular combustor with a design pressure loss ratio of 0. 055 at the design 
point. 




Specific fuel consumption, (kg/hr)/N 



Specific fuel consumption, (Ibm/h r)/lbf 


Figure 1. - Specific thrust plotted against specific fuel consumption for drone 
engine at cruise condition. Mach number, 0.8; altitude, 6100 meters 
(20 000 ft); pressure ratio difference, 0.99; compressor efficiency, 0.85; 
combustor efficiency, 0.96; turbine efficiency, 0.88; velocity coefficient, 
0.99; combustor AP/P, 0.96. 
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TABLE I. - STEADY-STATE PERFORMANCE 


(a) SI units 


Condition 


Altitude, m 
Mach number, M 
Ambient temperature, K 
Inlet recovery 
Speed, rpm 

Percent corrected speed 
Compressor 
Airflow, kg/sec 
Inlet temperature, K 

Outlet temperature, K 
2 

Inlet pressure, N/m 
Pressure ratio 
Adiabatic efficiency 
Bleed air, percent 
Combustor 

Airflow, kg/sec 

Outlet temperature, K 
2 

Outlet pressure, N/m 
Fuel flow, kg/hr 
Turbine 

Gas flow, kg/sec 

Outlet temperature, K 
2 

Outlet pressure, N/m 
Pressure ratio 
Adiabatic efficiency 
Nozzle 

Gas flow, kg/sec 
Inlet temperature, K 
Thrust coefficient 
Performance 
Thrust, N 

Specific fuel consumption, (kg/hr)/N 


10. 13x10 


38. 09x10 26. 54X10 


7. 09x10 
3. 96 
0.864 


1078. 3 


300. 78 


20. 68x10 
1.842 
0.85 


206.47 


14. 44x10 
1.84 
0.85 


2864. 5 
0. 107 


1547.9 
0. 133 


15. 29x10 15. 29x10 


1030.6 1098.9 663.3 

46. 68X10 4 49.62X10 4 14.37x10^ 
334.4 374.35 68.08 


25. 61x10 27. 22x10 


2224.0 
0. 150 


2664.8 
0. 141 











TABLE I. - Concluded. STEADY-STATE PERFORMANCE 


(b) U. S. Customary units 


Condition 

1 

2 

3 

4 

5 

Altitude, ft 

0 

20 000 

0 

0 

0 

Mach number, M 


0.8 

0. 8 

0.8 

0 

Ambient temperature, °R 

519.0 

447. 7 

519. 0 

519.0 

519.0 

Inlet recovery 

1.0 

0.99 

0.99 

0.99 

10 

Speed, rpm 

35 260 

34 800 

34 256 

35 200 

17 630 

Percent corrected speed 

100 

100 

91. 5 

94.0 

50.0 

Compressor 






Airflow, lbm/sec 

10.06 

7. 14 

12.86 

13.29 

4. 32 

Inlet temperature, °R 

519.0 

505. 1 

585. 1 


519.0 

Outlet temperature, °R 

810. 1 

786.6 

857. 2 


592.8 

Inlet pressure, psfa 

2116'. 8 

1481.8 

3193.9 


2116.8 

Pressure ratio 

3. 96 

3. 96 

3. 25 


1. 49 

Adiabatic efficiency 

0.859 

0.864 

0.861 

0.858 

0.849 

Bleed air, percent 

0.0 

0.0 

0.0 

0.0 

0.0 

Combustor 






Airflow, lbm/sec 

10.06 

7. 14 

12.86 

13. 29 

4. 32 

Outlet temperature, °R 

2010.0 

1941.0 

1855.0 

1978.0 

1194.0 

Outlet pressure, psfa 

7956.0 

5544. 0 

9748.8 

10 363.7 

3001.0 

Fuel flow, lbm/hr 

663. 1 

455. 2 

737. 2 

825.3 

150. 1 

Turbine 






Gas flow, lbm/sec 

10. 25 

7. 26 

13.06 

13. 51 

4. 36 

Outlet temperature, °R 

1752.0 

1701. 5 

1650.0 

1744.0 

1130.0 

Outlet pressure, psfa 

4318.6 

3015. 4 

5349.6 

5685. 1 

2345.8 

Pressure ratio 

1.842 

1.840 

1.822 

1.823 

1. 279 

Adiabatic efficiency 

0. 85 

0. 85 

0.85 

0.85 

0.85 

Nozzle 






Gas flow, lbm/sec 

10. 25 

7. 26 

13.06 

13. 51 

4. 36 

Inlet temperature, °R 

1752.0 

1701. 5 

1650.0 

1744.0 

1130.0 

Thrust coefficient 

0. 99 

0. 99 

0.99 

0. 99 

0. 99 

Performance 






Thrust, Ibf 

644.0 

mm 

500.0 


85.0 

Specific fuel consumption, (lbm/hr)lbf 


m 

1. 47 


1.77 









Engine System Simulation 


A mathematical description of the turbojet engine system dynamic and steady-state 
characteristics was formulated and used as the basis for an analog computer simulation. 
The engine simulation was then used for the investigation of various control schemes 
over a number of engine operating conditions. Since the system design was to be based 
on results obtained from the computer model, a comprehensive engine system simula- 
tion was developed. 

The overall engine system simulation can be regarded as a number of independent, 
but interconnected, simulations as shown in figure 2. External inputs to the engine sys- 


w f 



Figure 2. - Simulation diagram for ordnance engine. 


tem consist of compressor face conditions and the pressure into which the engine system 
exhausts. For this study, inlet dynamics were neglected and the compressor face total 
pressure and temperature were supplied as input constants for each operating condition. 
The details of the analog computer model are presented in appendix B. In addition to 
the analog model, a digital program was written in Fortran IV to verify the steady-state 
portion of the analog results. The digital model also provides a ready means for estab- 
lishing the effects of engine system modifications on steady-state performance. The 
digital program is presented in detail in appendix C. 


SYSTEM AND CONTROL DISCUSSION 

The computer simulation was first applied to determine the engine system steady- 
state operating characteristics. Once the steady-state performance was established, 
the model was used to investigate a number of possible fuel control systems. 
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Steady- State Characteristics 


The first simulation runs were used to define the normal engine operating character- 
istics. Figure 3 presents a compressor map obtained from the simulation. The com- 
pressor map was generated by holding rotor speed constant and varying the compressor 
exit impedance. The resulting variation in pressure ratio was plotted against corrected 
weight flow, and a constant speed operating line was thus obtained. The process was re- 
peated for a series of rotor speeds to form the complete map. As indicated in refer- 
ences 1 and 2, the compressor simulation technique applied in this system simulation 
can be used to define the compressor stall line. The stability boundary of the compres- 
sor simulation will closely match the compressor stall line. This stability boundary is 
included in figure 3. 

The engine normal operating line for 6100 meters (20 000 ft) and M = 0. 8 was also 
included in figure 3. This line is the locus of steady-state operating points plotted in 
the coordinates of the compressor map. The operating line shown in figure 3 illustrates 
engine performance at design conditions. As inlet conditions and exhaust nozzle area 
vary, the normal operating line will shift. 



2 3 4 5 6.7 8 9 10 11 

Corrected airflow, W 2 |/§ 2 / 62 , Ibm/sec 

Figure 3. - Compressor map for drone engine. 
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Engine normal operating lines were determined for a number of operating conditions, 
and are plotted on a steady- state compressor map in figure 4. The operating line for 
sea-level static conditions intersects the compressor stall line at 41. 5 percent corrected 
speed. The effect of inlet ram pressure, however, is to shift the operating line away 
from the surge line, such that there is sufficient surge margin at all Mach numbers 
greater than 0. 5 and corrected speeds greater than the windmilling speed. It must be 
recognized that the component data used in the computer model is preliminary. Final 
data for the engine model could shift the surge and normal operating lines. 



1 2 3 4 5 6 7 8 9 10 11 

Corrected airflow, Ibm/sec 

Figured - Compressor map with normal operating lines for sea-level and altitude conditions. 


The effect of varying the exhaust nozzle area at sea-level static conditions was in- 
vestigated. These results, shown in figure 5(a) indicate that additional surge margin, 
particularly at high rotor speeds, can be gained. However, the design point compressor 
performance varies with exhaust nozzle area. Figure 5(b) presents lines.of constant 
turbine inlet temperature superimposed on the compressor map. These data were ob- 
tained by varying fuel flow and exhaust nozzle area to establish an engine operating speed 
and a turbine inlet temperature. The normal operating line for the cruise condition is 
included for reference. 
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Stall line 



(a) Variation of normal operating line with exhaust nozzle area. 


Corrected turbine 
inlet temperature, 



altitude; Mach 0.8) 


i-i i i i i l I l 

0 1.5 2.0 2.5 3.0 3.5 40 45 5.0 

Corrected airflow, VI kg/sec 

i i i i i 1 I l l 

3 4 5 6 7 8 9 10 11 

Corrected airflow, Ibm/sec 

(b) Constant turbine inlet temperature lines. 

Figure 5. * Compressor map for drone engine. 







Net jet thrust, Ibf Actual airflow. Ibn/sec 



Corrected speed, N/y^, percent 

tcj Net jet thrust. (d) Specific fuel consumption 

Figure 7. - Drone engine performance at sea-level static condition. 
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Corrected turbine inlet temperature, T t 4/^ 2 Corrected net jet thrust, T^, Ibf 







The predicted propulsion system performance at Mach 0. 8 and 6100 meters 
(20 000 ft) is shown in figure 6. At 100 percent corrected speed, a turbine inlet temper- 
ature of 1078 K (1941° R), pressure ratio of 3. 96, and airflow of 3. 24 kilogram per sec- 
ond (7. 14 lbm/sec), the thrust (with a converging nozzle) is equal to 1548 newtons 
(348 lbf) with a specific fuel consumption of 0. 133 kilogram per hour per newton (1. 30 

(lbm/hr)/lbf). All performance predictions are for a nozzle exit area of 0. 0182 square 
2 

meter (0. 1963 ft ) and an inlet total pressure recovery of 1. 0. 

The system performance at sea-level static conditions is shown in figure 7. At 
100 percent corrected speed, a pressure ratio of 3. 96, and a turbine inlet temperature 
of 1117 K (2010° R), the thrust is 2865 newtons (644 lbf) with a specific fuel consumption 
of 0. 107 kilogram per hour per newton (1. 046 (lbm/hr)/lbf). 

Figure 8 presents the significant engine corrected performance parameters as a 
function of corrected speed for several flight conditions. These characteristics were ob- 
tained for an inlet total pressure recovery of 1. 0 and an exhaust nozzle area of 0. 0182 
square meter (0. 1963 ft^). - 


Windmilling Characteristics 

The simulation of the turbojet engine provided a convenient means of studying the 
free windmilling condition under various flight conditions. Of particular interest are the 
rotor speed, air flow, and combustor pressure and temperature. The windmill condition 
is achieved by simply reducing the fuel flow to zero and permitting the compressor and 
turbine torque to come to an equilibrium balance. 

It should be noted, however, that windmill conditions are difficult to estimate ac- 
curately. At windmill, both the compressor and turbine performance must be estimated 
far from the design point; in addition, the computer voltage levels are low, all of which 
introduce computational inaccuracies. 

For an engine which is to be started or restarted in flight from the windmilling con- 
dition, it is important to know the surge margin, air flow, and combustor pressure and 
temperature from the standpoint of combustor and igniter design and fuel scheduling. 

The combustion and igniter must be designed to ensure satisfactory ignition at the most 
severe windmilling condition anticipated, and the fuel-air mixture must be limited to 
values which preclude compressor stalling upon ignition. Further, the structural integ- 
rity of the rotor, particularly bearings, may be influenced by the windmilling condition. 

It was discovered that for all conditions investigated, the corrected windmill rotor 
speed is proportional to flight speed (fig. 9(a)). Corrected air flow through the engine, 
however, is a function only of flight Mach number (fig. 9(b)). Combustor pressure and 
compressor discharge temperature at the windmill condition are presented in figures 
9(c) and (d), respectively, for the anticipated range of flight conditions. 
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3050 <10 000) 
4570 (15 000) 
6100 (20 000) 
7620 (25 000) 




Hq ‘ajn^jaduja^aOJeyDSipjossaJduJOo 




I I I L/J 

aas/ujqi 'wijjib papaJJOD 


Flight Mach number Flight Mach number 

(b) Corrected airflow (c) Combustor pressure W) Compressor discharge temperature, 

plotted against plotted against Mach 

Mach number. number. 

Figure 9. - Windmill data for drone engine. 


Fuel Controls 


In addition to the investigation of basic engine performance characteristics, the sim- 
ulation was used to evaluate various engine fuel schedules and controls. Engine control 
concepts were examined for a specific application in which the engine would be started 
from the windmill condition in flight and, once started, would operate continuously at or 
near full thrust. 

The application requires a fuel schedule which supplies a prescribed amount of fuel 
as a function of the engine operating state. The fuel supplied must maintain engine oper- 
ation within safe operating limits during acceleration. The application does not neces- 
sarily require closed loop control and, if possible, should be satisfied by a simple fuel 
schedule. Several fuel controls were investigated with the system simulation. Three of 
these controls consisted of open loop fuel flow schedules, based on different engine pa- 
rameters. The fourth control was a generalized parameter fuel control, which was a 
more conventional closed loop system. 

The basic criteria used to evaluate the concepts was the ability to provide adequate 
acceleration within the compressor stall limit and to maintain the engine near design 
speed for cruise conditions between M = 0. 9 at 4570 meters (15 000 ft) and M = 0. 7 at 
7620 meters (25 000 ft). 

The selected fuel schedules are based on linear relations of a fuel flow parameter 
and compressor pressure ratio. These relations were arbitrarily selected using previ- 
ous control experience. 

Case A fuel control . - The case A fuel schedule was based on the corrected param- 
eter plots of figure 10. The curves plotted in this figure represent the system normal 
operating lines for various flight Mach numbers. These lines are represented in terms 
of a fuel flow parameter Wj/N 62 and the compressor pressure ratio P^. g/Pj. 2 * A 
point on a curve represents an equilibrium operating condition at a given flight Mach 
number, while a point above a normal operating line represents more fuel than required 
for equilibrium and, hence, system acceleration capability. 

The various normal operating lines converge to a single line as operation approaches 
the design compressor pressure ratio of 3. 96. It is evident from figure 10 that a linear 
schedule of the fuel flow parameter (dashed line) intersecting the normal operating line 
at the design point can be used to control the engine. This type of linear fuel schedule 
requires only the measurement of compressor inlet and discharge pressures and engine 
speed. The acceleration schedule must be above the normal operating lines to provide 
excess fuel to accelerate the system. The slope of the acceleration schedule was se- 
lected such that the engine is capable of accelerating from about 50 percent speed at 
sea-level conditions. 
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0 1 2 3 4 5 

Pressure ratio, P t 3 / P t 2 


Figure 10. - Engine corrected fuel flow to corrected speed requirements 
for various flight Mach numbers. 


Using these constraints, the fuel schedule was based on the relation 



These equations can be written as 

W f = N^O. 1311X10" 4 P t 2 + 0. 0414X10" 4 P t 3 ) kg/hr 

W f = N (l. 384X10" 3 P t 2 + °. 437X10" 3 P t 3 ) lbm/hr 

The aforementioned fuel schedule will accelerate the engine from windmilling to design 
speed and provide the correct fuel flow to maintain the engine at 100 percent corrected 
speed for all flight conditions. 
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0 12 3 4 


Pressure ratio, Pj 3 / P t 2 

Figure 11. - Engine corrected fuel flow requirements for 
various flight Mach numbers. 


Case B-fuel control. - The effects of scheduling the corrected fuel flow parameter 
W ( /6 2 V«2 as a function of compressor pressure ratio was also investigated. Figure 11 
presents normal operating lines for various flight Mach numbers plotted in terms of 
these parameters. It can be seen in figure 11 that the normal operating lines also con- 
verge to a single line as the design pressure ratio is approached. As with the case A 
fuel control, a linear fuel schedule through the design point and above the family of oper- 


ating lines can be used. 

The case B fuel control was then based on a linear schedule of V^2 against 

P. „/ P t 2 * The schedule for the case B control was adjusted to provide the correct fuel 
flow at 100 percent corrected speed. The equation for this control schedule is 


Wf 

P t 3 

= 2 4 + 74 91 - 

5 2^2 

P t,2 

t-h 

1 

P t 3 

5 29 + 165 14 ° 

6 2>/^2 

Is 0 Ll T JL.\J XJ rn li i 1 

P t,2 


kg/hr 


lbm/hr 
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or 


Wf = y[d ~ 2 (q. 2369X10 -4 P t 2 + 7 - 391X10' 4 P t 3 ) 

W f = y[f 2 ( 2 . 50x10 ” 3 P t 2 + 78. OlxlO" 3 P t 3 ) 

A direct mechanization of the fuel schedules given by equations (3) and (4) requires 
measurement of the compressor inlet total temperature. For the computer analysis, 
however, the j/ 0^ term was neglected. The schedule thus reduces to: 

W f = ( 0 . 2369X10" 4 P t 2 + 7. 391X10" 4 P t 3 ) kg/hr 

W f = (2. 50X10" 3 P t 2 + 78. OlxlO' 3 P t 3 ) lbm/hr 

The effects of neglecting the variation in compressor inlet temperature were examined 
during the computer study. The acceleration characteristics obtained with this fuel 
schedule are acceptable. The correct amount of fuel is supplied to maintain the engine 
within ±2. 5 percent of design speed of all flight conditions. The fuel schedule requires 
measurement of the compressor inlet and discharge pressure conditions. 

Case C fuel control . - A minor variation of the case B fuel schedule was also inves- 
tigated. From figure 11, it can be noted that the fuel flow intercept is extremely low 
and thus could be neglected. A linear fuel schedule of the following equation form was 
considered: 





or, again neglecting the temperature effects, 

W f = 0. 7350X10' 3 P, 

W f = 77. 57X10' 3 P t 

Generalized parameter control. - The controls discussed thus far schedule fuel flow 
as a function of engine parameters. The schedules presented result in engine accelera- 
tion and operation close to the design speed. Operation at other engine speeds, however, 
is not possible. The generalized parameter fuel control of reference 4 allows a more 
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complete engine control. This control was adapted for use with the drone engine, and 
simulated on the analog computer with the drone simulation. The generalized parameter 
fuel control consists of a series -parallel combination of orifices used to generate the de- 
sired control limits. The concept is based on the assumption that the maximum and mini- 
mum fuel limits can be described as linear functions of corrected engine parameters. 

The acceleration limit is defined by the following relation: 


W f N / P t 3 ^ 

= -^L- (0. 7439 — ^ + 0. 6804 




^2 V ^2 ' 


t, 2 


W 


l N 4. 64^,1. 


6 2\/^2 V^2 \ P t, 2 


50 


(kg/hr)/percent speed 


(lbm/hr)/percent speed 


> ( 10 ) 


and the minimum fuel limit is defined by 


W. 


— == = — | 0. 404 — ^ - 0. 404 + — 94x10 , 


6 2 V^2 V^2 


t,2 


t, 2 


4= = -4L= (o . 89 - 0. 89 + 1355 


6 2\^2 V ^2 


t,2 


t, 2 


(kg/hr)/percent speed 


MU) 


(lbm/hr)/percent speed 


J 


Dynamic Response of Engine and Fuel Controls 

A computer investigation of critical engine parameters during acceleration from 
windmill to full speed with the various fuel controls was made. The objective of this 
study was to establish the dynamic response of the engine system during the acceleration 
transient and to compare the different controls. 

The acceleration characteristics obtained from the engine and the generalized pa- 
rameter fuel control are presented first. As indicated in reference 4, this control con- 
cept results in engine dynamic performance similar to that obtained with conventional 
fuel controls. The data with this control is thus presented first to provide a basis for 
the evaluation of the other simplified fuel controls investigated. 

The operating characteristic obtained with the drone engine and the generalized pa- 
rameter control, during acceleration from windmill to design speed at the cruise condi- 
tion is presented in figure 12(a). For the large throttle advances represented by this 
transient, the governor spool travel is limited and acceleration occurs along the maxi- 
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mum fuel limit to the desired speed. The engine deceleration characteristic from design 
to 75 percent speed is also shown in figure 12(a). 

Similar acceleration and deceleration characteristics at sea-level static are pre- 
sented in figure 12(b). The acceleration characteristics for cruise and sea-level static 
are shown on a fuel parameter to pressure ratio plot in figure 13. The linear relations 
for the maximum and minimum fuel limits are evident in this figure. 



Figure 13. - Acceleration characteristics of drone engine with gener- 
alized parameter control. 


In general, the acceleration characteristics exhibited by each of the three pressure- 
related fuel controls (A, B, and C) were found to be similar. Figure 14(a) shows the 
acceleration lines for each of the three, and for the generalized parameter fuel control, 
plotted on a compressor map. The engine normal operating line and the compressor 
stall line are also included in the plot, thus giving a relative indication of the available 
acceleration and stall margins. The characteristics are plotted for the cruise condition 


21 




3 4 5 6 7 8 9 10 11 12 

Corrected airflow, Ibm/sec 

(b) Sea-level condition. 

Figure 14. - Acceleration characteristics for various control schedules. 
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of M = 0. 8 at 6100 meters (20 000 ft). The fuel schedules for the three simplified con- 
trols are nearly equivalent and the time to accelerate the engine from windmill to full 
speed was approximately 4 seconds for each of the three. The acceleration line for the 
generalized parameter fuel control crosses the lines generated by the simplified controls 
at about 70 percent corrected speed. The more conventional control (generalized param- 
eter) thus provides higher fuel flow above 70 percent speed and yields a somewhat faster 
acceleration transient. 

The acceleration characteristics for the four controls are presented for sea-level 
static conditions in figure 14(b). At this operating condition the engine normal operating 
line and the acceleration lines are closer to the compressor stall line in the low speed 
range than they were at the cruise condition. Because of the resulting low surge mar- 
gin at the windmill condition, engine accelerations were made from an idle condition at 
60 percent corrected speed. The general trends, however, were found to be similar to 
those obtained at the cruise condition. 

The transient acceleration performance of the drone engine was investigated at the 
cruise altitude condition. As an example, the fuel schedule B was selected to illustrate 
the transient performance from windmilling to design speed. Typical responses of sig- 
nificant engine variables are presented in figure 15. Immediately after ignition, the 




(a) Speed, fuel flow, and actual airflow transients. 

Figure 15. - Acceleration characteristics with fuel schedule B at 6100-meter (20 000-ft) altitude and Mach 0. 8 from 
windmilling to design speed. 
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(b) Pressure transients. 



0 12 3 

Time, sec 

(c) Temperature transients. 


Figure 15. - Concluded. 



compressor discharge pressure and temperature increase. The turbine inlet tempera- 
ture reaches 820 K (1475° R) and then rises to the final value of 1078. 3 K (1941° R) as 
the engine accelerates to 100 percent corrected speed. 


Effect of Cruise Altitude on Fuel Control Performance 

The simulation results were used to predict the effects of various fuel controls on 
flight dynamics of a missile having a nominal drag force of 1557 newtons (350 Ibf) at a 
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Figure 16. - Engine, control, and vehicle flight dynamics. 
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Mach number of 0. 8 and an altitude of 6100 meters (20 000 ft). The engine and fuel con- 
trol are designed to produce 1557 newtons (350 lbf) thrust at this flight condition. 

Figure 16 shows the variation of thrust and drag with Mach number for altitudes of 
4570, 6100, and 7620 meters (15 000, 20 000, and 25 000 ft, respectively). The thrust 
plotted is that which would result from the use of the fuel control calibrated to give 
1557 newtons (350 lbf) nominal thrust at the nominal cruise condition of Mach 0. 8 at 
6100 meters (20 000 ft) altitude. Drag is assumed to be typical of a missile designed to 
cruise at this flight condition. 


CONCLUSIONS 

The computer analysis of the drone engine indicates that the engine performance is 
satisfactory over the desired flight range. It is shown that sufficient stall margin is 
available to accelerate the engine from windmilling to design speed. 

The fuel schedules A, B, and C will accelerate the engine over the specified range 
of flight conditions to design speed. The generalized parameter control provides satis- 
factory performance over the desired flight range of the engine. This type of control is 
recommended for the drone application. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, January 13, 1972, 

132-15. 
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APPENDIX A 


SYMBOLS 

A area, m 2 ; ft 2 

b coefficient, kg/hr; lbm/hr 

coefficient 

C n coefficients, (rpm) 2 /(K) 2 ; (rpm) 2 /(°R) 2 

Cp specific heat at constant pressure, J/(kg)(K); Btu/(lbm)(°R) 

E energy, N-m/sec; ft-lbf/sec 

g gravitational constant, 1 (kg)(m)/(N)(sec 2 ); 32. 17 (lbm)(ft)/(lbf)(sec 2 ) 
h enthalpy, J/kg; Btu/lbm 

Ah enthalpy drop, J/kg; Btu/lbm 

2 2 

I inertia of rotating mass, N-m-sec ; lbf-ft-sec 

J mechanical equivalent of heat, 1 N-m/J; 778.3 ft-lbf/Btu 

K n coefficient, (rpm)(sec)/m; (rpm)(sec)/ft 

k coefficient, (N 2 )(sec 2 )/(kg 2 )(m 4 )(K); (lbf 2 )(sec 2 )/(lbm 2 )(ft 4 )(°R) 

k^ coefficient, (N)(sec)/(kg)(K) 4//2 (m) 2 ; (lbf)(sec)/(lbm)(°R) */ 2 (ft) 2 

Lq compressor torque, N-m; ft-lbf 

L t turbine torque, N-m; ft-lbf 

l length, m; ft 

M Mach number 

N rotational speed, rpm 

P pressure, N/m 2 ; lbf/ft 2 

R universal gas constant, 287 (N)(m)/(kg)(K); 53. 3 (lbf)(ft)/(lbm)(°R) 
r mean radius, m; ft 

rrp tip radius, m; ft 

s Laplace operator, sec - * 

T temperature, K; °R 

t time, sec 
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U m mean rotor speed, m/sec; ft/sec 

V volume, m^; ft^ 

v„ axial velocity, m/sec; ft/sec 

z 

W weight flow, kg/sec; lbm/sec 

Wj* externally acted upon fuel flow, kg/sec; lbm/sec 
a coefficient, 30 sec/min 

/3 rotor air inlet angle, deg 

y ratio of specific heats, 1. 4 

6 ratio of total pressure to standard atmospheric pressure 

7 j efficiency 

9 ratio of total temperature to sea-level atmospheric temperature 

3 3 

p weight density, kg/m ; lbm/ft 

(p flow coefficient 

p 

\f/ pressure coefficient 

T 

temperature coefficient 
Subscripts: 
a actual 

c variable associated with stage characteristics 

d design conditions 

f fuel 

1 ideal 

n stage number designation 

t total condition 

tr total condition reference state 

tv total condition variable associated with stage volume 

v variable associated with the stage volume 

0 free-stream condition 

2 compressor inlet 

3 compressor discharge 
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APPENDIX B 


ENGINE SYSTEM SIMULATION - ANALOG REPRESENTATION 

The analog computer representation of the drone engine is based on the engine simu- 
lation techniques presented in references 1 and 2. The basic simulation consists of 
representations of the engine components and their associated gas dynamics. The com- 
ponent representations are interconnected to form the complete engine system model. 


Compressor Simulation 


The modeling techniques used to represent the four-stage axial compressor are 
based primarily on the method developed in reference 1. Each stage is modeled with 
steady- state stage characteristics and the gas dynamics of an associated stage volume. 
The steady- state performance of the individual compressor stages is represented by 
pressure and temperature rise coefficients i// and »£/ _ which are functions of a flow 

C ' C 

coefficient <p. These coefficients are mathematically defined for the n tn stage as 


* 


P 

c, n 


T 

n 1 tv, n- 1 
N 


P \ 2 / 7 
tc,n \ 


tv, n-1/ 



(Bl) 


IT 


tv, n-1 


(B2) 


<P 


c, n 


s 

N 


zc, n 


(B3) 


where 



P T 

Figure 17 presents plots of the pressure and temperature coefficients \p Q and \{/ c as 
functions of the flow coefficient cp, for the four stages of the engine of this report. 

The flow coefficient cp is proportional to the ratio of the axial flow velocity and the 
rotor speed N. The axial-flow velocity can be related to corrected weight flow by 



The angle /3, defined as the angle between the air velocity vector and the axial direction, 
was assumed small for the engine system considered. In the actual computer model, the 
complexities of equation (B4) were avoided by the use of a quadratic fit to a numerical 
evaluation of the equation. The steady-state representation of the compressor was thus 
based on the pressure and temperature rise coefficients of figure 17 and the approxima- 
tion of equation (B4). 
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Flow function, 

(b) Temperature rise coefficients. 

Figure 17. - Compressor stage characteristics. 
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The gas dynamics for each compressor stage were assumed to occur in an equivalent 
stage volume conceptually located just downstream of the stage pressure and temperature 
rise characteristics. The stage dynamics were represented by the following continuity, 
energy, and momentum balances across the equivalent stage volume: 


1 <Psv ; n)=^-<W c ,n-W C;n+1 ) 


dt 


d 


dt ( p sv, n T tv, n^ y ( T tc, n W c, n " T tv, n W c, n+1^ 
v n 


_d_ 

dt 


(W „ Js^(P, 


tc, n 


- P tv,n) 


(B5) 

(B6) 

(B7) 


To compute the total pressure in the equivalent stage volume, an equation of state of the 
following form was assumed: 


tv, n 


P ^sv, n P tv, n^ 


(B8) 


The stage characteristics, the continuity, energy and momentum balance equations, 
and the equation of state form a complete set for a single compressor stage. Figure 18 
presents a block diagram illustrating the computational sequence used for a typical 
stage. The stage inlet pressure, temperature, and flow are combined to determine the 
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c,n 


p tv,n-l 

p tv,n-l 



p tv,n 
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Figure 18. - Block diagram of simulation for compressor stage. 
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flow coefficient cp. The flow coefficient is then used to determine the stage pressure 
and temperature rise coefficients, and thus the discharge total pressure and tempera- 
ture. The discharge conditions from the stage maps and the lumped volume gas dy- 
namics are used to compute the conditions in the stage volume, and hence the inlet con- 
ditions for the succeeding stage. The overall representation for the complete four-stage 
compressor was formed by interconnecting the individual stage models. 


Combustor Simulation 

The combustor simulation was based on a lumped parameter representation of the 
system gas dynamics and a steady-state map for the efficiency of the combustion proc- 
ess. The combustion efficiency was assumed to be the nonlinear function of pressure 
and temperature rise presented in figure 19. The simulation did not include combustor 
blowout limits. 

Momentum, continuity, and energy balances were used to represent the system 
lumped parameter gas dynamics. A detailed discussion of the combustor equations and 
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Combustor parameter, P t 4 AT ( 4, (lbf/ft^)°R 

Figure 19. - Combustor efficiency representation. 
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simulation techniques is presented in reference 2. The pressure drop across the burner 
was determined from 


, W^T t 3 
AP' = k - 
’ P 

*t,3 


(B9) 


and assumed to be 5 percent of the compressor discharge pressure at the design point. 

The combustion process itself includes some delays. The combustion delays were 
included in the simulation by arbitrarily adding a delay of the following form between 
commanded fuel flow and the fuel flow supplied to the combustor: 


W” e -0. 005 s 
W { (1 + 0. 003 s) 


(BIO) 


Turbine Simulation 

The steady-state representation of the single-stage turbine was based on the curves 
of figures 20 and 21. The efficiency curve of the first figure was mechanized directly 



Figure 20. - Turbine efficiency. 
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Figure 21. - Turbine flow characteristics. 


on a computer function generator, while the more complex flow to pressure ratio rela- 
tion of the second figure was approximated by 


where 



(Bll) 


The turbine efficiency characteristic is given as a function of rotor speed and actual 
enthalpy drop across the turbine. Functionally the relation can be described by 


^ 5 



(B12) 
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where 


In the computational sequence, the actual enthalpy drop is assumed and then used to 
determine the turbine efficiency from the map of equation (B12). The actual enthalpy 
is also used to determine the turbine discharge temperature from 

< T t, 4 - T t, 5> = jr 4h a ( B13 > 

The discharge temperature is then used in a continuity balance, across the volume just 



T t, 4 


Figure 22. - Computational diagram for tu rbine. 
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downstream of the turbine, to compute the discharge pressure. In equation form, this 
is 


RT 


t.5 


t, 5 


/ < w 5 - 


(B14) 


The discharge pressure is in turn combined with turbine inlet conditions to compute the 
ideal enthalpy drop from 


Ah i " C p T t, 4 I 



(y-l)/y-| 


(B15) 


The ideal enthalpy drop is multiplied by the efficiency to determine the actual drop as- 
sumed earlier in the computational sequence. The calculations are shown in schematic 
form in figure 22. 


Exhaust Nozzle Simulation 

The engine exhaust nozzle, located just downstream of the turbine exhaust, was the 
final engine element included in the simulation. Flow through the fixed area nozzle was 
computed from 



Rotor Dynamics Simulation 


The engine rotational speed is determined by integrating a torque balance across 
the compressor and turbine 
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(B17) 



(B18) 


(B19) 


(B20) 


and the engine speed equation then becomes 

— - = — Vl q + Lrp) (B21) 

dt Ul / C 1 

where I is the mass polar moment of inertia of the rotor assembly. 

The various component representations were interconnected to form the engine sys- 
tem simulation. The system simulation was then used to investigate and define engine 
operating characteristics. 
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APPENDIX C 


ENGINE SYSTEM SIMULATION - DIGITAL REPRESENTATION 

The digital program uses the system equations to calculate steady-state solutions 
at various percentages of design speed. 

The simulation includes: the four-stage compressor, the burner, the turbine, and 
the nozzle. The MAIN program calls subroutines which may be classified into three 
groups: the input-output subroutines, the system subroutines, and the data handling 
subroutines. 

Input-Output Subroutines 

INPUT sets constants used in other subroutines 
OUTPUT prints desired variables 

System Subroutines 

CMPRSR solves compressor equations for each stage 
BURNER calculates burner variables 
TURBIN calculates turbine variables 
NOZZLE calculates nozzle variables 

Data-Handling Subroutines 

TERPl two-dimensional interpolation for table lookups 

TERP2 three-dimensional interpolation for table lookups 

BLDATA block data subroutine, includes all necessary data for table lookup 

INPUT sets in the area constant for the nozzle, ambient pressure, temperature, 
specific heat of compressor and burner, gamma of compressor and burner, and various 
physical constants. Mach number is also calculated but never used. 

CMPRSR is entered with weight flow (estimated), pressure, temperature, speed, 
and geometric and physical constants. It calculates the outlet pressure and temperature 
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(using phi, psi, and psi prime maps), efficiency, change in enthalpy, and outlet weight 
flow for each stage. The final calculations are the overall efficiency of the compressor 
and the sum of the enthalpy changes of the four stages of the compressor. 

BURNER is entered with the weight flow, pressure and temperature of the compres- 
sor exit, the estimated outlet temperature, and various constants. The enthalpy and the 
outlet pressure are calculated. The product of the outlet pressure and the delta temper- 
ature is calculated and is used in a table lookup to find efficiency. Fuel flow is then 
calculated. 

TURBIN is entered with the weight flow, pressure and temperature (estimated) of 
the burner exit, compressor speed, fuel flow, the torque of the compressor, and vari- 
ous constants. A table lookup is used to find the efficiency, and from that the outlet 
pressure and temperature are calculated. 

NOZZLE is entered with weight flow, pressure and temperature of the turbine exit, 
and various constants. Nozzle weight flow, outlet pressure and temperature, and noz- 
zle exit weight flow are calculated. 

OUTPUT calculates compressor pressure ratio, corrected weight flow, and cor- 
rected speed. It then prints all variables of interest. 

TERP1 is a function subprogram which performs two-dimensional interpolation. 

TERP2 is a function subprogram which performs three-dimensional interpolation. 

BLDATA is a block data subroutine which contains data for table lookups of psi, 
psi prime, phi, tip radii, mean radii, areas, P t ^ AT, burner efficiency, N/j/T t 
P t 4/Pt 5, weight flow parameter, turbine enthalpy parameter, and turbine efficiency. 

A flow chart for the drone engine is presented in figure 23. The Newton-Raphson 
method of convergence is used in both the burner temperature iteration and the inlet 
weight flow iteration. 
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Figure 23. - Flowchart for drone engine. 
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Digital Computer Program 


1 1 B FTC MAIN LIST 

COMMON / VAR / N, PI, TI, WI, SUMDHC, ETACOM, THET, T(4), P(4), 

1 W ( 4 ) , ETAC ( 4 ) , 0ELHCI4), PHI ( 4 ) , PS I ( 4 ) » P S I P t 4 ) , MC ( 4 ) , WFG » 

2 DHT , ETAT, P3, T3, W3, T4, P4, W4, T5, P5, W 5, T6, P6, W6 , H7, 

3 ET A B , T4G 

COMMON / CONST / GAMMA, GAMMA I , GAMMAB, G, CPC, CPB, AJ, RK, R2P » 
1 SCCNST, PA, TA, MA, HC , C0NST7 

COMMON / BLOA / XPHI(30,4), YPSI(30,4), YPSIP(30,4), R(4I, RM(4), 

1 A { 4 ) , SPEEOII2), XP40TBI 16 ) , YET AB ( 16 ) , XP4P5(19), XNT4(3), 

2 YW5RTPI 19,31, XUTERMI 25) , YEFFI25) 

REAL N, MC, MA, NT4 


PERCENT OF RADII t 100.0 
GEOMR = 1.0 

C SQUARE ROOT OF PERCENT OF AREA / 100.0 

GEOMA = 1.0 

C PRINT PAGE CONSTANT 

L = 0 

WRITE (6,3) 

3 FORMAT (1H1): 

DO 2 J = 1 ,4 
R ( J ) = GEOMR * R ( J ) 

RM ( J ) = GEOMR * RM ( J ) 

2 A ( J ) = GEOMA * GEOMA * A(J) 

C INLET FLOW, PRESSURE, AND TEMPERATURE ESTIMATES 

WI = 15.0 
PI = 20.4 
TI = 518.7 

C BURNER TEMPERATURE ESTIMATE 

T4G = 1960.0 
CO 30 K = 1,10 
WI = .93 * WI 
CALL INPUT 

THET = SORT (T I / 518.7) 

XX = THET * SPEED! K ) / 100.0 
N = 35260.0 * XX 

SCONST = 19.7392 *N*N/AJ/G/ CPC / 3600.0 / 144.0 
IJCCNV = 0 
1 CALL OUTPUT ( L ) 

CALL CMPRSR 
KCONV = 0 
LGCCNV = 0 

5 CALL BURNER 
T 4 = T4G 

6 CALL TURBIN 

7 CALL NOZZLE 

C CONVERGENCE OF WEIGHT FLOWS IN TURBINE AND NOZZLE 

ERW = W5 - W7 

IF (ABS(ERW) .LE. .0001 * W5) GO TO 15 
LGCCNV = LGCONV + 1 
IF (LGCONV .LT. 50) GO TO 36 
WI = .999 * WI 
GO TO 1 
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36 IF (ERW .NE. ERWP) GO TO 35 
WI = .999 * WI 

GD TO 1 

35 IF (KCONV .GT. 0) GO TO 10 
TAGF = TAG 
ERWF = ERW 
TAG * .99 * TAG 
KC 0T> V = 1 
GO TO 5 

10 TAGN = ( TAG * ERWP - TAGP * ERWI / (ERWP - ERW) 

IF (TAGN .GT. T3) GO TO 37 

TAG = 1.25 * T3 
WI = 1.001 * WI 
GO TO 1 

37 IF ( TAGN .NE. TAG) GO TO 3A 
WI = .999 * WI 

GO TO 1 

3A IF (TAGN .LT. .99 * TAG) TAGN = .99 * TAG 

IF (TAGN .GT. 1.01 * TAG) TAGN = 1.01 * TAG 

TAGP = TAG 

ERWF = ERW 

TAG = TAGN 

KCONV = KCONV ♦ 1 

GO TO 5 

C CONVERGENCE OF INLET AND OUTLET WEIGHT FLOWS 

15 NT A = N / SQRT(TA) 

IF ( NT A .GT. XNTA(l)) GO TO A 
IXN = 2 
XXN = 0.0 
GO TO 11 

A CO 5 IXN = 2,3 
IXN = IXN 

IF (NTA .GE. XNTA(IXN-l) .AND. NT A .LT. XNTA(IXN)) GO TO 12 
9 CONTINUE 
IXN = 3 
XXN = 1.0 
GO TO 11 

12 XXN = (NTA - XNT A ( I NX- 1 ) ) / (XNTA(IXN) - XNTA(IXN-l)) 

11 PAP5 = PA / P5 

IF (PAP5 .GT. XPAP5(1)) GO TO 13 
IXP = 2 
XXP = 0.0 
GO TO 17 

13 CO 1A IXP = 2,19 
IXP = IXP 

IF ( PAP5 .GE. XPAP5 ( IXP-1 ) .AND. PAP 5 .LT. XPAP5UXP)) GO TO 16 
1A CONTINUE 
IXP = 19 
XXP = 1.0 
GO TO 17 

16 XXP = ( PAP 5 - XPAP5( IXP-1 ) ) / ( XPAP5 ( IXP) - XP AP5 ( I XP- 1 ) ) 

17 W5RTP = TERP2IXXP, XXN, YW5RTP ( I XP , I XN ) , YW5RTP ( I XP-1 , I XN) , 

1 YW5RTP( IXP, IXN-1) , YW5RTP( IXP-1, IXN-1) ) 

W5C = (W5RTP * PA / SQRT(TA)) 

WOUT = W5C “ WFG 

IF (AeSIWI - WOUT) .LE. .0001 * ABS(WI)) GO TO 25 
IF (IJCONV .GT. 0) GO TO 20 
WIP = WI 
WOUT P = WOUT 
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WI = 1.005 * WI 
IJCCNV = 1 
GO TO 1 

20 WIN = ( WOUTP * WI - WOUT * WIP) / ((WI - WIP) - (WOUT - WOUTP)) 


IF 

(WIN 

• IT. 

.994 

* WI) 

WIN = .994 * 

WI 

IF 

(WIN 

• GT. 

1.006 

* WI ) 

WIN = 1.006 

* WI 

IF 

(WIN 

• GT. 

WI ) 

T4G = 

(WI / WIN) * 

T4G 

IF 

(WIN 

• LT. 

WI ) 

T4G * 

(WIN / WI ) * 

T4G 

WIP 

= WI 






WOUTP = 

WOUT 





WI 

= WIN 







IJCCNV = IJCONV * 1 
GO TO 1 

25 CALL OUTPUT ( L ) 

30 CONTINUE 
RETURN 
END 


$ I BFTC INPUT LIST 

SUBROUTINE INPUT 

COMMON / CONST / GAMMA* GAMMA 1 * GAMMAB, G, CPC, CPB, AJ, RK, R2P, 
1 SCCNST, PA, TA, MA, HC , C0NST7 
REAL MA 


C0NST7 = 56.3 

MA * .8 

PA = 1 A . 7 

TA = 518.7 

CPC = .24 

CPB * .278 

GAMMA = 1.4 

GAMMAB = 1.34 

G = 32.16 

AJ = 778.12 

RK = 640.32 

GAMMA1 = GAMMA - 1.0 

R2P = RK f 6.2832 

HC = 18640.0 

RETURN 

end 
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$ I BFTC CMPRSR LIST 

SUBROUTINE CMPRSR 

COP NON / VAR / N, PI, TI, WI, SUMDHC, ETACOM, THET, T ( A ) , P(A), 

1 W(A), ETAC(A), DELHC(A), PHI ( A ) , PSIIA), PSIP(A), MC(A), WFG , 

2 DHT, ETAT, P3, T3, W3, TA, PA, WA, T5, P5, W5 , T6, P6, W6 , W7, 

3 ETAB, TAG 

COP P ON / CONST / GAMP1A, GAMMA 1 , GAMMAB, G, CPC, CPB, AJ, RK, R2P 
1 SCCNST , PA, TA, MA, HC , CONST7 

COPPON / BLDA / XPHI ( 30, A ) , YPSI(30,A), YPSIP(30,A), R ( A) , RM ( A ) 

1 A ( A ) , SPEED! 12) , XPADTB (16), YETAB ( 16 ) , XPAP5(19), XNTA<3), 

2 YW5RTP( 19,3) , XUTERM(25), YEFF ( 25 ) 

REAL N, MC, MA 


STAGE 1 

WC = WI * SQRTIRK * TI) / PI / A(l) 

PCONV = 0 

VG = WI * RK * TI / PI / A( 1) 

1 RAD = VG * VG / RK / T I * GAMMA 1 / 2A.0 / G / GAMMA 
IF (RAD .LE. 1.0) GO TO 6 

WI = 100.0 * WI 
RETLRN 

6 WTERM = VG / SORT ( RK * TI) * ((1.0 - RAD) ** (1.0 / GAMMA1 ) ) 
ERW = WC - WTERM 

IF (ABS(ERW) ,LE. .000001 * ABS(WC)) GO TO 3 

IF (MCONV .GT. 0) GO TO 2 

VGP = VG 

ERWF = ERW 

VG = .99 * VG 

PCONV = 1 

GO TO 1 

2 VGN = (VG * ERWP - VGP * ERW) / ( ERWP - ERW) 

VGP = VG 

ERWP = ERW 
VG = VGN 

PCONV = MCONV + 1 
GO TO 1 

3 PHIII) = 60.0 * VG / 6.2832 / RM(1) / N 
IF (PHI(l) .GT. XPHI(1,H) GO TO A 

IX = 2 
GO TO 10 

A CO 5 IX = 2,30 
IX = IX 

IF IPHI(l) .GE. XPHI(IX-1,1) .AND. PHI(l) .LT. XPHI(IX,1)) 

1 GO TO 10 
5 CONTINUE 
IX = 30 

10 XX = (PHIII) - XPHI ( IX-1, 1 ) ) / ( XPH I ( I X , 1 ) - XPHI ( IX-1,1) ) 

15 PS I ( 1 ) = TERPMXX, YPSKIX.l), YPSI ( I X-l ,1 ) ) 

IF ( PS I ( 1 ) .NE. 0.0) GO TO 103 
CO 101 L = 1,30 
K = 32 - L 

IF (YPSI (K-1,1) .NE. 0.0) GO TO 102 

101 CONTINUE 

102 XK = (PHI(l) - XPHI(K-1,1) ) / ( XPH I ( I X , 1 ) - XPHI(K-1,1)) 

PS I ( 1 ) = TERP1IXK, YPSI(K,1), YPSI(K-1,1)) 

103 PSIF(l) = TERPKXX, YPSIP(IX,1), YPS I P ( I X- 1 , 1 ) ) 

IF (PSIP(l) .NE. 0.0) GO TO 106 



DO 1 04 L = 1,30 
K = 32 - L 

IF (YPSIPIK-1,1) .NE. 0,0) GO TO 105 

104 CONTINUE 

105 XK = (PHI(l) - XPHI (K-l , 1 ) ) / ( XPHI ( I X , 1 ) - XPHI(K-l,i)) 

PSIP(l) = TERPKXK, YP S IP ( K » 1 ) , YPSIP ( K-l, 1 ) ) 

106 ETAC(l) = PSIP(l) / PSI(l) 

MC(1) = SQRT ( ( 9. 8696 * RM(1) * RM(1) * N * N / TI - VG * VG / TI) 

1 / (21600.0 * G * RK * (1.0 - GAMMA1 * VG * VG / 24.0 / G / RK / 

2 GAMMA / TI ) ) ) 

P ( 1 ) * P I * ((PSIP(l) * SCONST * R ( 1 ) * R ( 1 ) / TI ♦ 1.0) ** 

1 (GAMMA / GAMMA1)) 

T ( 1 ) = Ti + PS I ( 1 ) * SCONST * R ( 1 ) * R(l> 

DElFC(l) = WI * CPC * ( T ( 1 ) - TI) 

W( 1 ) = WI 

C STAGES 2 THROUGH 4 

CO 35 I = 2,4 
W ( I ) * W(I-l) 

WC = W(I) * SQRT ( RK * T(I-l)) / P(I-l) / A U) 

NCONV = 0 

VG = W(I) * RK * T(I-l) / PCI-1) / All) 

16 RAD = VG * VG / RK / T(I-l) * G4MMA1 t 24.0 / G / GAMMA 
IF (RAD .LE. 1.0) GO TO 36 

WI = 100.0 * WI 
RETURN 

36 WTERM = VG / SQRT( RK * TCI-l)) * ((1.0 - RAO) ** (1.0 / GAMMAD) 
ERW = WC - WTERM 

IF (ABS(ERW) .LE. .000001 * ABS(WC)) GO TO 18 

IF ( NCONV .GT. 0) GO TO 17 

VGP = VG 

ERWF = ERW 

VG = .99 * VG 

NCONV = 1 

GO TO 16 

17 VGN = (VG * ERWP - VGP * ERW) / (ERWP - ERW) 

VGP = VG 

ERWP = ERW 
VG = VGN 

NCONV = NCONV + 1 
GO TO 16 

18 PH I ( I ) = 60.0 * VG / 6.2832 / RM(I) / N 
IF IPHIU) .GT. XPH 1(1,1)) GO TO 19 

IX = 2 
GO TO 25 

19 DO 20 IX = 2,30 
IX = IX 

IF ( PH I ( I ) .GE. XPHI ( IX — 1, 1 ) .AND. PHI(I) .LT. XPHI(IX,D) 

1 GO TC 25 

20 CONTINUE 
IX = 30 

25 XX = (PHI(I) - XPHI ( IX-1, I ) ) / ( XPHI ( I X , I ) - XPHI ( I X-l , I ) ) 

30 PSIU) * TERPKXX, YPSI(IX,I), YPS I ( I X-l , I ) ) 

IF (PSI(I) .NE. 0.0) GO TO 109 
DO 107 L = 1,30 
K = 32 - L 

IF (YPS I (K-l , I ) .NE. 0.0) GO TO 108 

107 CONTINUE 

108 XK = ( PH I ( I ) - XPH I ( K-l , I ) ) / ( XPH I ( I X , I ) - XPHI(K-1,I)) 

PSIU) = TERPKXK, YPSI(K,I), YPSI(K-1,I)) 
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109 PSIPtn = TERP1(XX, YPSlPUX.I), YPSIPtIX-1, 1 )) 

IF (PSIP(I) .NE. 0.0) GO TO 32 

CO 110 L = 1,30 
K s 32 - L 

IF ( YPS I P( K-l , I ) .NE. 0.0) GO TO 111 

110 CONTINUE 

111 XK = ( PH I ( 1 ) - XPHI(K-l,m / ( XPHI ( I X , I ) - XPHI(K-1,D) 

PSIF(I) = TERPUXK, YPSIP(K,I ), YPSIP(K-1,I)) 

32 ETACI I ) - PSIP(I) / PSI(I) 

33 *C(I> = SORT ( (9.8696 * RM ( I ) * RM(I) * N * N / T(I-l) - V G * VG / 

1 T(I-l)) / (21600.0 * G * RK * (1.0 - GAMMA 1 * VG * VG / 24.0 / G 

2 / PK / GAMMA ! T( I-D) )) 

P(I) = P(I-l) * ((PSIP(I) * SCONST * R(I) * R ( I ) / T(I-l) + 1.0) 

1 ** (GAMMA / GAMMAD) 

T ( I ) = T(I-l) ♦ PSI(I) * SCONST * R ( I ) * R(I) 

CEOCU) = W(I) * CPC * ( T ( I ) - TU-ll) 

35 CONTINUE 
P3 = P(4) 

T3 = T(4> 

W3 = W ( 4 ) 

SUMCHC =0.0 
CO 40 I = 1,4 

40 SUMCHC = SUMDHC + DELHC(I) 

ETACCM = ((P3 t PI) ** ( GAMMA 1 / GAMMA) - 1.0) / <(T3 / TI ) - 1.0) 

RETURN 

END 
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tIBFTC BURNER LIST 

SUBROUTINE BURNER 

CORRON / VAR / N, PI, Tit WI , SUMDHC , ETACOM, THET, T(4), P(4), 

1 W(4), ETAC (4 ) , DELHC ( 4 ) , PHI ( 4 ) , PSII4), PSIP(4), MC(4), WFG, 

2 DHT, ETAT, P3, T3, W3, T4, P4, W4, T5, P5, W5, T6, P6, W6 , W7, 

3 ETAB, T4G 

CO PR ON / CONST / GAMMA, GAMMA I , GAMMAB, G, CPC, CPB, AJ, RK, R2P, 
1 S_CC NS T , PA, TA, MA, HC , C0NST7 

CORRON / BLDA / XPHI(30,4), YPSI(30,4), YPSIP(30,4), R(4), RM(4), 

1 A ( 4 ) , SPEED! 12 ) , XP4DTB(16), YET AB ( 16 ) , XP4P5(19), XNT4I3), 

2 YW5RTP! 19,3) , XUTERM ( 25 ) , YEFF ( 25 ) 

REAL N, MC, MA 


W4 = W3 

1 P4 = P3 - .00214 * W4 * W4 * T3 / P3 
P4DTB = P4 * I T4G - T3) ♦ 1800.0 

IF ( P4DTB .GT. XP4DTB(1)) GO TO 2 
IX = 2 
XX = 0.0 
GO TO 15 

2 CO 5 IX = 2,16 
IX = IX 

IF ( P4DTB .GE. XP4DTBUX-1) .AND. P4DTB .LT. XP4DTBUX)) GO TO 10 
5 CONTINUE 
IX * 16 
XX = 1.0 
GO TO 15 

10 XX = (P4DTB - XP4DTBI IX-1 ) ) / (XP4DTBUX) - XP4DTB ( I X-l ) ) 

15 ETAE = TERPKXX, YETAB(IX), YETAB(IX-l)) * .9915 
H4 = T4G * CPB 

WFG = W4 * (CPB * T3 - H4) / (H4 - ETAB * HC ) 

RETURN 

END 
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$ I 8FTC TUP B I N LIST 

SUBROUTINE TURBIN 

COUPON / VAR / N, PI, Tit WI« SUMDHC , ETACOM, THET, T(4), P(4), 

1 W C 4 ) , E TAC ( 4 ) , DELHC ( 4 ) , PHK4I, PSI ( 4 ) , PSIP(4), MC(4), WFG » 

2 DHT, ETAT, P3, T3, W3, T4, P4, W4, T5, P5, W5, T6, P6, W6, W7, 

3 ETAB, T4G 

COMRON / CONST / GAMMA, GAMMA 1 , GAMMAB, G, CPC, CPB, AJ, RK, R2P 
I SCCNST, PA, TA, MA, HC , CONST7 

CORN ON / BLDA / XPHI(30,4), YPSI(30,4), YPSIP(30,4), R(4), RM(4) 

1 A ( 4 ) » SPEED ( 12 ) , XP4DT B ( 1 6 ) , YET AB ( 16 ) , XP4P5(19), XNT4(3), 

2 YW5 RTP (19,3), XUT ERM ( 2 5 ) , YEFF ( 25 ) 

REAL N, MC, MA, NT4 


W5 = W4 + WFG 
CHACT = SUMCHC / W 5 
DHT = SUMDHC 
U = 8.73E-3 * 3.97 * N 
UT = U * U / G / AJ 
UTEPM = UT / DHACT 

IF (UTERM .GT. XUTERMI1)) GO TO 50 
IXU = 2 
XXU = 0.0 
GO TO 65 

50 CO 55 IXU = 2,25 
IXU = IXU 

IF (UTERM .GE. XUTERM ( I XU-1 ) .AND. UTERM .LT. XUTERM (IXU)) 

1 GO TO 60 
55 CONTINUE 
IXU = 25 
XXU = 1.0 
GO TO 65 

6C XXU = (UTERM - XUT ERM ( I XU- 1 ) ) / (XUTERM(IXU) - XUTERM ( I XU- 1 ) ) 

65 ETAT = TERPKXXU, YEFF(IXU), YEFF(IXU-l)) 

CHIC = DHACT / ETAT 

P5 = P4 * ((1.0 - DHID / CPB / T4G ) ** (GAMMAB / (GAMMAB - 1.0))) 

T5 = T4G - DHACT / CPB 

RETURN 

END 



$ I BFTC NOZZLE LIST 

SUBROUTINE NOZZLE 

COPPON / VAR / N, PI, TI, WI, SUMDHC , ETACOM, THET, T<4), P(4), 

1 W(4), ETAC ( 4 ) , OELHC ( 4 ) , PHI(4), PSI ( 4) , PSIP(4), MC ( 4 ) , WFG , 

2 DHT, ETAT, P3, T3, W3, T4, P4, W4, T5, P5, W5 , T6, P6, W6 , W7, 

3 ETAB, T4G 

COP PON / CONST / GAMMA , GAMMA1 , GAMMAB, G, CPC, CPB, AJ, RK, R2P, 
1 SCCNST , PA, TA, MA, HC , C0NST7 
REAL N, MC, MA 
C 

C 

W6 = W5 
P6 = P5 
T6 * T5 

P A P( = PA / P6 

IF (PAP6 .LT. .5283) P AP6 = .5283 
IF (PAP6 .GT. 1.0) PAP 6 = 1.0 

V»7 = C0NST7 * P6 * (PAP6 ** (1.0 / GAMMA)) * SORT (1.0 - ( PAP6 ** 

1 (GAPPA1 / GAMMA))) / SORT ( T6 ) 

RETLRN 

END 
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$ I BFTC OUTPUT LIST 

SUBROUTINE OUTPUT(L) 

COVNON / VAR / N, PI, TI, WI, SUMDHC, ETACOM, THET, T(4), P(4), 

1 W(4), ETAC ( 4 ) , OELHC ( 4 ) , PHI ( 4 ) , PSII4), PSIP(4), MC(4), WFG, 

2 DHT, ETAT, P3 f T3, W3, T4 f P4, W4, T5, P5, W5, T6, P6 f W6 , W7 f 

3 ET AB , T4G 
REAL N, MC 


L = L + 1 
POPI = P3 / PI 

COR V - WI * THET / (PI / 14.7) 

CORN = (N / THET / 35260.0) * 100.0 

WRITE (6,10) N, WFG, PI, TI, WI, P(l), T(l), W(l), ETAC(l), 

2 OELHC(l), P ( 2 ) , T ( 2 ) , W(2), ETAC(2), DELHC ( 2 ) , P(3), T<31, W(3), 

3 ET AC ( 3 ) , OELHC ( 3 ) ♦ P(4), T(4), W(4), ETAC(4), DELHC ( 4 ) , P3, T3, 

4 W3, ETACOM, SUMDHC, P4, T4, W4, ETAB, P5, T5, W5, ETAT, DHT, P6, 

5 T6 , W6 , W7 , POPI, CORW , CORN 

10 FORMAT (1H0* 9HN = , 1PE14.6, 10X, 9HWFG = , 1PE14.6 / 

1 IX, 7HSTAT ION , 16X, 8HPRES SURE , 17X, 4HTEMP , 21X, 1HW, 22X, 

2 3HEFF, 21 X , 2H0H / IX, 5HINLET, 4X, 1P3E24.6 / IX, 6HC0MP 1, 3X, 

3 1P5E24.6 / IX, 6HC0MP 2, 3X, 1P5E24.6 / IX, 6HC0MP 3, 3X, 

4 1P5E24.6 / IX, 6HC0MP 4, 3X, 1P5E24.6 / IX, 9HC0MP EXIT, 1P5E24.6 

5 / IX, 6HBURNER, 3X, 1P4E24.6 / IX, 7HTURBINE, 2X, 1P5E24.6 / IX, 

6 6HN0ZZLE, 3X, 1P3E24.6 / IX, 8HN0Z EXIT, 49X, 1PE24.6 / IX, 

7 7HF0PI = , 1PE14.6, 10X, 11HC0RR W * f 1PE14.6, 10X, 12HPERCENT 

8 N = , 1PE14.6) 

WRITE (6,20) ( PHI ( I ) , PSI(I), PSIP(I), MC ( I ) , I * 1,4) 

20 FORMAT ( 2 3 X , 3HPHI, 21X, 3HPSI, 21X, 4HPSIP, 15X, 

1 1 3E INLET MACH NO / IX, 6HC0MP 1, 1P4E24.6 / IX, 6HC0MP 2, 

2 1PA624.6 / IX, 6HC0MP 3, 1P4E24..6 / IX, 6HC0MP 4, 1P4E24.6) 

IF (L .LT. 3) RETURN 

WRITE (6,30) 

30 FORMAT ( 1 H 1 ) 

L = 0 

RETURN 

END 


TERP1 LIST 

FUNCTION TERP1 (X, A, B) 


TERF1 =B+(A-B)*X 

RETURN 

END 
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$ I BF TC TEFP2 LIST 

FUNCTION TERP2 (X, Y, A, B, C, 0) 


TERF2 =D+<B-D)*Y+(C-D)*X-(C+B-A-D)*X*Y 

RETLRN 

END 


IBFTC BLCATA LIST 
BLOCK DATA 

COMMON / BLDA / XPHI(30,4), YPSI<30,4), YPSIP(30,4), R<4), RM<4), 

1 MM. SPEED! 12 ) * XP4DTBI 16 ) , YETABI 16) , XP4P5U9), XNT4I3), 

2 YW5RTP( 19,3) , XUTERM(25), YEFF ( 25 ) 


PSI DATA 

DATA ((YPSI(I,K)» I = 1,30), K = 1,4) / 

C PSI PRIME DATA 

DATA ( I YPS I P ( I , K ) , I = 1,30), K = 1,4) / 

C PHI DATA 

DATA < ( XPH I ( I »K ) , I = 1,30), K = 1,4) / 

C TIP RADII DATA 

CATA R / 

C MEAN RADII DATA 

DATA RM / 

C STACE AREAS DATA 

DATA A / 

C DESIRED PERCENTAGES OF DESIGN MECHANICAL SPEED 

DATA SPEED / 

C P4DELTA T PARAMETER DATA 

DATA XP4DTB / 

C BURNER EFFICIENCY DATA 

CATA YETAB / 

C N/T4 PARAMETER DATA 

DATA XNT4 / 

C P4/F5 PARAMETER DATA 

DATA XP4P5 / 

C TUREINE WEIGHT FLOW PARAMETER DATA 

DATA (<YW5RTP<I,J), I = 1,19), J = 1,3) / 

C TUREINE TORQUE PARAMETER DATA 

DATA XUTERM / 

C TURBINE EFFICIENCY DATA 

DATA YEFF / 

END 
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